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FEELD OF THE INVENTION 

This invention relates to numerical simulation, and more specifically relates to an 
improved method of simulating a hydrocarbon system that may include one or more 
10 of the following: a subterranean hydrocarbon-bearing formation, inj ection wells that 
penetrate the formation, production wells, surface flowlines, associated aquifers, and 
surface fluid processing facilities. 

BACKGROUND OF THE INVENTION 

Numerical simulation is widely used in industrial fields as a method of simulating a 
15 physical system by using a computer. In most cases, there is desire to model the transport 
processes occurring in the physical system. What is being transported is typically mass, 
energy, momentum, or some combination thereof. By using numerical simulation, it is 
possible to model and observe a physical phenomenon and to determine design 
parameters, without actual laboratory experiments and field tests. 
20 Reservoir simulation is of great interest because it infers the behavior of a real 

hydrocarbon-bearing reservoir from the performance of a model of that reservoir. The 
typical objective of reservoir simulation is to understand the complex chemical, physical, 
and fluid flow processes occurring in the reservoir sufficiently well to predict future 
behavior of the reservoir to maximize hydrocarbon recovery. Reservoir simulation often 
25 refers to the hydrodynamics of flow within a reservoir, but in a larger sense reservoir 
simulation can also refer to the total hydrocarbon system which can include not only the 
reservoir, but also injection wells, production wells, surface flowlines, associated 
aquifers, and surface processing facilities. Reservoir simulation calculations in such 
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hydrocarbon systems are based on fluid flow through the entire hydrocarbon system being 
simulated. These calculations are performed with varying degrees of rigor, depending on 
the requirements of the particular simulation study and the capabilities of the simulation 
software being used. 

5 The principle of numerical simulation is to numerically solve equations describing 

a physical phenomenon using a computer. Such equations are generally ordinary 
differential equations and partial differential equations. As a means for numerically 
solving such equations, there are known the finite element method, the finite 
difference method, the finite volume method, and the like. Regardless of which 

10 method is used, the physical system to be modeled is divided into cells (a set of which 
is called a grid or mesh), and the state variables that vary in space throughout the 
model are represented by sets of values for each cell. The reservoir rock properties 
such as porosity and permeability are typically assumed to be constant inside a cell. 
Other variables such as fluid pressure and phase saturation are specified at specified 

15 points, sometimes called nodes, within the cell. A link between two nodes is called a 
"connection." 

Since the reservoir simulation can include vastly different fluid flow environments 
(e.g., porous rock, well tubing, processing facilities), the set of cells can include multiple 
segments of different flow enviroments. Although individual segments, such as 
20 production facilities and surface pipe segments, could be represented by single cells, at 
times reservoir simulation programs subdivide such segments into multiple cells. 

A set of equations is developed to express the fundamental principles of 
conservation of mass, energy, and/or momentum within each cell and of movement of 
mass, energy, and/or momentum between cells. These equations can number in the 
25 millions. The replacement of state variables that vary in space throughout a model by 
a finite number of variable values for each cell is called "discretization". In order to 
analyze a phenomenon changing in time, it is necessary to calculate physical 
quantities at discrete intervals of time called timesteps, irrespective of the 
continuously changing conditions as a function of time. Time-dependent modeling of 
30 the transport processes therefore proceeds in a sequence of timesteps. The simulation 
of flow during a desired period of time is done by incrementally moving forward 
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through time, and solving the equations at each timestep for the pressure and fluid 
contents of every node. 

Efforts have been made to perform widely varying reservoir simulation methods 
in a single computer code. However, such "general-purpose" simulation systems are 
5 very complex, largely because they are designed to have one or more of the following 
capabilities: 

(1) represent many different types of cells (e.g., different domains of a reservoir, 
well tubing, and surface gathering and distribution facilities, and surface 
processing facilities); 

10 (2) use different timestep computation methods (e.g., IMPES, fully implicit, 

sequential implicit, adaptive implicit, and/or cascade methods); 

(3) use different ways of representing reservoir fluids; 

(4) use multiple ways of computing transport between cells; 

(5) perform what is called a "black-oil model," which treats the hydrocarbons as 
15 being made up of two components, and also having the capability of 

performing compositional representations in which the hydrocarbons are 
assumed to contain compounds such as methane, ethane, propane, and heavier 
hydrocarbons; 

(6) simulate steam injection or in situ combustion processes, which must take into 
20 account temperature changes as a function of time which require an energy 

balance and related calculations; 

(7) simulate miscible recovery processes using special fluid properties and 
transport calculations; 

(8) simulate hydrocarbon recovery processes that take into account injection of 
25 surfactants, polymers, or other chemicals and the flow of these fluids into the 

reservoir; 

(9) simulate injection of chemicals that react with each other, the reservoir 
hydrocarbons, or the reservoir rock; and 

(10) simulate migration of hydrocarbons and geologic deposition over geologic 
30 time. 

Historically, because general-purpose simulation systems are complex, reservoir 
simulation programs typically either model only fluid movement in the subsurface 
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reservoir or they model the subsurface fluid movement and include simplified 
calculations, for example, hydraulics models, to model flow within well tubulars 
between the reservoir and the surface. To more realistically model the flow through 
the wells and in attached surface production facilities (such as manifolds, pumps, 
5 compressors, separators and pipelines) has typically required a separate facilities 
simulation model that is distinct from the reservoir simulation model. Since these 
reservoir and facilities simulation models attempt to describe the flow behavior over a 
long period of time, typically many years, and the simulations are discretized into 
timesteps that cover the desired simulation period, it is cumbersome and costly to try 

10 to devise an integrated modeling process based on two distinct models. Such a 

process usually involves alternating timesteps between the two models, such that the 
changes that occur in one model can be incorporated into the input for the subsequent 
timestep of the second model. 

One challenging problem in reservoir simulation in the past has been the lack of 

15 rigorous modeling of transport phenomena through well tubulars and surface facilities 
within a reservoir simulation system. Conventional reservoir simulators typically 
have two significant limitations. First, many reservoir management studies require a 
more sophisticated surface facilities flow model than is possible in a reservoir 
simulator. To address this problem, the simulator user typically develops and 

20 maintains one simulation model for the reservoir behavior and another simulation 

model for the facilities behavior, with no convenient way to integrate the results from 
the two models, or to easily evaluate the effects of change in one model on the 
behavior of the other. Second, reservoir management typically involves dramatic 
changes in well or facility rates. Since facility rates are used as boundary conditions 

25 for reservoir simulation calculations, dramatic changes in those boundary conditions 
often lead to numerical instabilities in the simulation calculations that are difficult to 
manage. 

Another challenging problem in reservoir simulation is the ability of a simulator 
user to specify facility management logic in the simulation. To optimize the 
30 production of hydrocarbons from a field requires a continuous process to monitor and 
control the production and injection facilities during the life of the field, which can 
span up to 30 or more years. To more accurately predict the behavior of the reservoir 
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using a reservoir simulator requires that the impact of surface facility controls be 
taken into consideration during the simulation. The specific controls, and their 
sequence and timing, is difficult, if not practically impossible, to program into the 
simulator because the permutations would be limitless and depend on each field and 
5 each simulator user's approach to managing the field. Therefore, the reservoir 

simulator user needs a tool to be able to specify to the simulator the custom checks 
and controls that the simulator user believes will optimize the production of the 
reservoir. Those checks and controls then need to become an integral part of 
simulator computations. In short, an intuitive tool is needed to capture each simulator 

10 user's custom facility management logic, turn that into executable code that is 

integrated with the reservoir simulator's code, and then execute the resulting system to 
predict the overall behavior of the reservoir during its lifetime. 

The typical approach to this problem has been the development of software that 
provides a user with a structured capability. The simulator user is provided with a 

15 text-based editor or interface wherein he or she has available a set of constructs and a 
fixed set of keywords. Using the constructs, the simulator user can define potential 
conditions, which if met during the simulation, trigger some response that modifies 
the behavior of the wells. This approach is not very intuitive for the user. The 
process of assembling the constructs and keywords can be tedious, and the result is 

20 not easy to read or interpret, particularly for persons not originally involved with the 
reservoir project. There is no clear overview of the logic, meaning that it is not easy 
to follow the high-level line of reasoning behind the logic. The user has limited 
control over the process and the remedies available for execution are limited by a 
finite set of combinations that the simulator recognizes, resulting in limited flexibility 

25 of the system to accommodate unconventional reservoir management strategies. 

There is a need in the industry for a reservoir simulation system that (1) provides a 
more realistic model of flow through well tubulars and surface facilities and (2) 
captures a facility management logic customized by a reservoir simulator user and 
executes it as part of a reservoir simulation. 
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SUMMARY 

The invention relates to a computer system and method for simulating transport 
phenomena in a complex system. The computer system comprises a logic interface 
that enables a user of the computer system to dynamically construct logic to 

5 customize simulation of the physical system, a means for converting the constructed 
logic into corresponding object-oriented code, a means for integrating the object 
oriented code with the main simulation system which comprises a simulation data 
model and simulation algorithms, resulting in an integrated simulation system, and a 
means for executing the integrated simulation system using a plurality of connected 

10 processors. 

In carrying out the method of the invention, the computer system user carries 
out the steps of dynamically constructing logic to customize simulation of the 
physical system, initiating simulation of the physical system causing initiation of the 
following steps (i) automatically converting the logic into corresponding object- 

15 oriented code, (ii) integrating the object oriented code with the main simulation 

system which comprises a simulation data model and simulation algorithms, resulting 
in an integrated simulation system for simulating the physical system, and (iii) 
executing the integrated simulation system using a plurality of connected processors 
to simulate the physical system. 

20 The invention is particularly useful in simulating a reservoir system. In a 

reservoir system model, a facility network model ("FNM") is preferably developed 
that extends the discretized reservoir simulation model (consisting of nodes and flow 
connections) beyond the reservoir to include nodes and connections for modeling 
fluid flow in the well tubulars and surface production and injection facilities (such as 

25 manifolds, pumps, compressors, gathering lines, separators and pipelines). Each 
distinct facility type is modeled as a specialized type of node, connection, or a 
composite of several nodes and connections. In its most basic form, the entire 
simulation model can encompass every component of the hydrocarbon field from the 
subsurface reservoir, all wells and well hardware, and surface facilities up to and 

30 including the product delivery outlet(s). 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The new reservoir simulation system and its advantages will be better understood 
by referring to the following detailed description and the attached drawings. 

Fig. 1 illustrates a discretized reservoir simulation model and its attached network 
5 of wells and surface facilities. 

Figs. 2 and 3 illustrate a generic C++ class hierarchy showing a Facility Network 
Model. Fig. 3 is part of the class hierarchy shown in Fig. 2. 

Fig. 4 illustrates a simple facility network as a graphical representation of icons 
and arrows labeled with the name of the facility that each represents. 
10 Fig. 5 illustrates a Logic Diagram that shows a simplified facility management 

logic. 

Fig. 6 illustrates a schematic of a small portion of an example simulator class 
hierarchy. 

Fig. 7, extending the diagram in Fig. 6, illustrates how dynamic specialization is 
15 used to support user-written operations. 

The drawings are not intended to exclude from the scope of the new reservoir 
simulation system other embodiments that are the result of normal and expected 
modifications of these specific embodiments. 

DETAILED DESCRIPTION 

20 A computer system, having memory means, storage means, and an object-oriented 

software product can be used in simulating complex systems such as two- and three- 
dimensional domains of the hydrocarbon system that are discretized into structured 
grids, unstructured grids, or a combination of both. The computer system can also be 
used in situations in which the computational approach yields a topology having more 

25 than three dimensions, such as occurs in simulating fluid flow through fractured 
porous media. The simulation system is useful in simulating a characteristic of a 
complex system in which a transport phenomenon is occurring. The term "transport 
phenomena" is used in this description in a broad sense to include momentum 
transport (viscous flow), energy transport (heat conduction, convection, and 

30 radiation), and mass transport (diffusion). For example, the simulation system is 
particularly useful to simulate removal of hydrocarbons from reservoir rock using 
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various supplemental recovery techniques such as, but not limited to, thermally-based 
recovery methods such as steam flooding operations, waterflooding operations, and 
gas-drive-based methods that can operate under either miscible or immiscible 
conditions. The simulation system can also be used in modeling contaminant 

5 transport through subterranean formations. In this description, the simulation system 
will be described with respect to performing a general-purpose reservoir simulation 
on unstructured grid systems, and is especially useful for parallel computing in 
simulating hydrocarbon-bearing systems, which can include not only the reservoir, 
but also associated aquifers, production wells, injection wells, surface flowlines, 

10 manifolds, separators, valves, pumps, compressors, or any other surface processing or 
transport facility. 

For the purposes of the detailed description of the present invention, several terms 
are defined to clarify their meaning. A "reservoir simulation model" is a specific 
mathematical representation of a real hydrocarbon reservoir. An engineer managing a 

15 hydrocarbon reservoir may create many different reservoir simulation models, 
possibly with varying degrees of complexity, in order to quantify the past 
performance of the reservoir and predict its future performance. In the present 
invention, a reservoir simulation model is a discretized representation of a reservoir's 
physical space and the wells and surface facilities that are used in the recovery of 

20 hydrocarbon fluids. A "facility" as used in this description is a representation of a 
tangible piece of physical equipment through which hydrocarbon fluids are either 
produced from a reservoir or injected into a reservoir. In its broadest sense, the term 
facility is applied to any equipment that may be present along the flow path between a 
reservoir and its delivery outlets, which are the locations at which hydrocarbon fluids 

25 either leave the model (produced fluids) or enter the model (injected fluids). Facilities 
may comprise production wells, injection wells, well tubulars, wellhead equipment, 
gathering lines, manifolds, pumps, compressors, separators, surface flow lines and 
delivery outlets. In some instances, the term "surface facility" is used to distinguish 
those facilities other than wells. A "facility network" is the complete collection of 

30 facilities that are present in the model, which would include all wells and the surface 
facilities between the wellheads and the delivery outlets. 



EXPRESS MAIL NO. ET607889812US 
DECEMBER 6, 2001 



9 



A reservoir simulation model comprising a three-dimensional hydrocarbon- 
bearing reservoir and its associated facilities is illustrated in Fig. I. The reservoir 
simulation model schematically illustrated in Fig. 1 is comprised of a discretized 
reservoir 100, which is penetrated by wells 106a, 108a and 1 10a, depicted by dashed 

5 lines, and containing additional facility nodes 102, 103 and 104 and facility 

connections 1 12, 1 14, 1 16, 1 18 and 120. In Fig. 1, the reservoir 100 is depicted as 
comprising a plurality of volumetric cells 101, delineated by narrow solid lines. 
Although not explicitly shown in Fig. 1, each cell 101 represents a three dimensional 
volume of space that contains one node at its center. Also not explicitly shown, each 

10 pair of neighboring reservoir nodes has a connection between the corresponding 
nodes. In the reservoir simulation model, hydrocarbon fluids flow through a 
connection between a pair of nodes based on the pressure difference between the 
nodes and the transport properties of the connection. Thus, hydrocarbon fluids can be 
transported from a volumetric cell to a neighboring volumetric cell. Each well 106a, 

15 108a and 1 10a contains one node each, 106b, 108b and 1 10b, respectively. Also not 
explicitly shown in Fig. 1, each well's node is attached via a connection to each of the 
reservoir nodes through which the corresponding well's trajectory passes. The facility 
nodes 103 and 104 and facility connections 112, 114, 116, 118 and 120 represent a 
gathering network of flowlines that carry fluids produced from reservoir 100 into 

20 wells 106a, 108a and 1 10a to a common delivery point 102. 
Obiect-Oriented Programming 

Implementation of the present invention uses object-oriented programming 
("OOP") techniques. The most common programming languages for OOP are 
Simula, Eiffel, C++, Smalltalk, Objective-C, or variations thereof. However, the 

25 practice of the present invention is not limited to a particular OOP language. C++ is 
used in this detailed description to help the reader better understand and appreciate the 
benefits of OOP. The design of OOP is well known to those skilled in the art of OOP 
systems and will only be described generally. 

In an object-oriented program, the focus is primarily on data and secondarily on 

30 the functions that access that data, rather than being primarily on functions and 
secondarily on the data they require. In contrast with a program that is described 
mainly by procedures with data input to and output from each procedure (such as a 
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program written in FORTRAN or COBOL), an object-oriented program is organized 
around "objects." An object is a data structure and a set of operations or functions 
that can access that data structure. OOP systems will typically contain a large number 
of objects. Each operation (function) that can access the data structure is called a 
5 "method." The data that an object contains are called "member variables" and the 
methods that an object has are called "member functions." An object can be 
considered a smart collection of data, because the object's member functions give the 
object the ability to make changes to itself, that is, to change the values of its member 
variables. 

10 A "class" is a template for a specific type of object. The class defines the 

information associated with the object (its data) and the operations it can perform (its 
methods). When an object is created, the class is used as a template to define what 
data and methods the object has. An object is said to be an "instance" of the class that 
defines it. There can be many distinct objects, all instances of the same class, but 

15 each instance has its own data and methods. 

A typical object-oriented software design contains at its core a group of classes, 
which may be referred to collectively as a "data model." A given class in a data 
model defines one specific type of "object," which can contain data ("member 
variables") and methods ("member functions"). The data model is a major part of the 

20 object-oriented software's "architecture." An object-oriented software product's 
architecture encompasses the complete design of the software, including the data 
model around which it is built. 

OOP enables objects to model virtually any real-world entity in terms of its 
characteristics, represented by data, and its behavior, represented by the methods it 

25 can perform using the data. In this way, objects can model concrete things like cells 
in a reservoir simulation and abstract concepts like numbers. For the purposes of the 
present invention, benefits of OOP arise out of three basic principles: encapsulation, 
inheritance, and polymorphism. 

Objects hide, or encapsulate, the internal structure of their data, i.e., the explicit 

30 collection of member variables that an object contains, and the algorithms by which 
their member functions work. Instead of exposing these implementation details, 
objects present interfaces that represent their abstracts clearly with no extraneous 
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information. The interface allows only a limited number of interactions between the 
outside world and an object. In C++ programming, most of these interactions involve 
invoking member functions of the object. By invoking an object's member functions, 
the outside world can tell the object to do something, but the details of how the object 

5 performs this action are not visible to the outside world. In effect, the object 
encapsulates its data and methods and hides them from the outside world. 

For example, an object-oriented data model for a software system that calculates 
node properties might have a class called "BaseNode." The BaseNode class could 
have member variables such as "x_location," "yjocation," and "z_location" which 

10 hold numeric values to define a BaseNode object's position in space, and member 
functions to change the value of the BaseNode object's x, y, or z coordinate value. A 
software system built upon this example data model might contain code to create and 
modify any number of distinct BaseNode objects. Each BaseNode object contains its 
own unique values for its member variables, and each has the built-in member 

15 functions to modify those member variables. From the perspective of the software 
system, it does not matter how a BaseNode object maintains its location information. 
The software system only needs BaseNode member functions to set and retrieve this 
data. This organization of data (member variables) and code (member functions) into 
classes is referred to as encapsulation. 

20 Inheritance allows software developers to design the data model in the form of a 

"class hierarchy," such that related classes mimic the natural classification of object 
types. The three main benefits of inheritance are that it permits software developers 
to design code that encapsulates data efficiently, it permits implementation of 
methods on a specific class so that these methods can be used for objects of many 

25 specialized types, and it promotes reuse of pre-existing design and code. In a class 
hierarchy, general data and methods belong to a "base" class or "superclass." More 
specialized classes ("subclasses") "inherit" the base class data and methods and define 
additional data and methods of their own. These specialized classes are "derived" 
from the base class. Inheritance allows a software developer to define classes and the 

30 objects which are later created from them and related through the class hierarchy. 
Specifically, classes may be designated as subclasses of other base classes. A 
subclass inherits and has access to all of the functions of its base classes as though 
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such functions appeared in the subclass. Inheritance allows previously written 
programs to be broadened by creating new superclasses and subclasses of objects. 
New objects are described by how they differ from preexisting objects so that entirely 
new programs need not be written to handle new types of data or functions. 
5 For an example of a class hierarchy that illustrates inheritance, extending the 

previous example of the BaseNode class, a second class called "WellNode" might 
inherit from the BaseNode class. A WellNode object is both a WellNode and a 
BaseNode, but the WellNode class can define data and code specific to a WellNode. 
If a third class "NetworkNode" is defined, that also inherits from BaseNode, but not 

10 from WellNode, then a NetworkNode object will not contain the WellNode-specific 
data or code. However, both a WellNode object and a NetworkNode object will 
contain BaseNode- specific data and code. This simple example only shows a 
BaseNode class with one level of additional specialization (the WellNode and 
NetworkNode classes). However, inheritance permits any number of specialization 

15 layers, such that a WellNode class may have further subclasses, and the BaseNode 
class may have further superclasses. 

Polymorphism takes encapsulation a step further. A software component can 
make a request of another component without knowing exactly what that component 
is. The component that receives the request interprets it and determines, according to 

20 its variables and data, how to execute the request. Using the simple class hierarchy 
described previously as an example, a software component can invoke a member 
function of a BaseNode object. The software component does not need to know 
whether the BaseNode object is actually a WellNode object or a NetworkNode object. 
However, the class hierarchy enables the WellNode object and the NetworkNode 

25 object to contain different implementations of the invoked member function. Each 
class in the hierarchy may contain the same abstract interface for the member 
function, but they can have different functionality, appropriate to the specific type of 
object that executes the member function. 
Facility Network Model 

30 The reservoir simulation system of the present invention preferably includes a 

Facility Network Model ("FNM"). The FNM is designed to minimize the need for 
code changes and recompilation. This is achieved by defining generic classes in the 



EXPRESS MAIL NO. ET607889812US 
DECEMBER 6, 2001 



13 



FNM C++ class hierarchy. A representative generic C++ class hierarchy is shown in 
Figs. 2 and 3. Each box in Fig. 2 (numbered 200 through 218) represents a C++ class. 
Solid lines between classes with an open triangle at one end indicate inheritance 
relationships between classes, with the base class at the end with the triangle and the 

5 derived class at the opposite end. Dashed lines between two classes represent 
associations between objects of those two classes. The text labeling a dashed line 
indicates the name of the association, and the numeric labels at the two ends of each 
association indicate the cardinality of the association. For example, the association 
"node_is_child_of ' 220 has a "1" label at the "Node" end and a M 0..*" label at the 

10 "Connection" end. These labels indicate that every Connection object has exactly 1 
Node that is that Connection object's child, and every Node object has zero-or-many 
Connection objects that the Node object is a child of. Each generic class contains 
only the minimal set of member variables and member functions to define generic 
behavior for all objects of that class. For example, a generic node class (Node 203 in 

15 Fig. 2) is the most specialized type of class in the class hierarchy for all facility types 
that are modeled by some specialization of a node, whether that node represents a 
location within a well, surface flow line, a separator, or some other facility that can be 
represented by a node. Separate generic classes are used to define named "attributes" 
for specialized facility types (see Attribute Value 21 1 and its derived classes 212 and 

20 213). Using generic facility type classes (200 through 205) and generic named 

attribute classes (207 through 213) allows essentially all of the specialized behavior of 
desired facility types to be defined in a "Definitions Data File" (referred to hereafter 
as "DDF"). Appendix 6, at the end of this description, sets forth an example of this 
file. The key design feature of using generic facility and generic attribute classes 

25 enables the FNM to be flexible and extensible so that new facility types or new 

attributes can be incorporated through additions to or modifications of the contents of 
the DDF. This feature avoids additions or changes to the classes in the data model, 
which would require significant modifications to source code and recompilation of 
modified code. 
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The FNM data model is continued in Fig. 3, which shows a series of classes that 
keep track of the facility types and attribute types that are defined in the DDF. 

The DDF contains text descriptions for each possible facility type, such as 
"WellNode," "NetworkNode, " "SeparatorNode," "Well," and "NetworkConnection." 
5 It also contains text descriptions for each of the specialized named attributes for those 
facility types. When the DDF is processed by the software, it creates a 
Facility TypeSystemAttributes object (303) for each specialized facility type in the 
DDF, and an AttributeDefmition object (304) for each named attribute in the DDF. A 
FacilityTypeSystemAttributes object inherits from base classes VarBase 300, 

10 Variable AttributesBase 301 and Facility Type 302. The facility has system attributes 
association (305) relates a specific FacilityTypeSystemAttributes object (the definition 
for a particular facility type) to all of the AttributeDefmition objects (the definition for 
a specific named attribute) for that facility type. After the DDF has been completely 
processed, the reservoir simulation system contains objects that define the possible 

15 facility types (the collection of FacilityTypeSystemAttributes objects 303) and 

attributes (the collection of AttributeDefmition objects 304) that each of those facility 
types must have. The definition of available facility types and named attributes 
remains constant in the software until the DDF is modified and processed again, which 
would typically occur only as new facility types or attributes are needed. 

20 The reservoir simulation system having a FNM contains a graphical user interface 

(GUI) for defining the network of facilities ("Facility Network") that are part of a 
given simulation model. The GUI displays the Facility Network as a collection of 
icons that are attached to one another via connection arrows. 

Fig. 4 illustrates this graphical user interface for a simple Facility Network. 

25 Different icons can be used to represent distinct facility types (wells, separators, 

network junctions), and the connection arrows represent flowlines between the "node- 
like" facilities. In Fig. 4, each text box (with "b" suffixes 401b through 410b) 
indicates the name of a facility depicted by icons or arrows (with "a" suffixes 401a 
through 41 0a). The example Facility Network in Fig. 4 contains three wells 406a, 

30 408a, 4 1 0a and the surface flow network to which these wells are attached. The 
surface network contains NetworkNodes 403a and 404a, TerminalNode 401a, 
NetworkConnections 405a, 407a and 409a, and TerminalConnection 402a. Wells 
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406a and 408a produce fluids that move to the NetworkNode 404a via the 
NetworkConnections 405a and 407a. Fluids from Well 410a merge with the 
production from the other two wells at NetworkNode 403 a. From NetworkNode 
403a, the fluids continue to the delivery outlet TerminalNode 401a. 

5 The FNM GUI provides custom icons, arrows, dialogs, menus, and toolbar 

buttons and text so that the user, using one or more of these features, can build-up, 
edit, and visualize a Facility Network with relative ease. When a user adds a specific 
facility to the Facility Network for a reservoir simulation study, the reservoir 
simulation system uses the appropriate Facility TypeSystemAttributes object 303 

10 (given the type of facility that the user is adding) to locate the set of 

AttributeDefmition objects (304) which provides the system with the information it 
needs to create the actual facility object (202, 203, 204 or 205) and all of its system- 
defined attribute objects (212). 
Support of Custom Facility Management Logic 

15 The generic C++ data model for the FNM facilitates the complementary custom 

Facility Management Logic system ("FML") described below. The FML system is 
designed for monitoring and controlling the operation of facilities, including 
mechanical facilities associated with the field, within a running simulation model to 
optimize the production of hydrocarbons during the life of a field. The FML operates 

20 on facility objects during a simulation run that are created from and analogous to the 
FNM data model objects. To maximize the customizability of the FML, the FNM 
data model incorporates user-defined named attributes (class UserAttributeValue, 
213, Fig. 2). A user can create additional named attributes for any of the available 
facility types (beyond those that are defined by users in the DDF) and then use those 

25 custom named attributes within their custom FML. The definition for a user-defined 
named attribute is saved as an AttributeDefmition object (218, Fig. 2), and is 
associated with a Facility TypeUser Attributes object (217, Fig. 2) that corresponds to 
the facility type for which the attribute has been defined. 

The FNM can be used with any suitable object-oriented programming language. 

30 The data model design that includes generic facility types, generic named attributes, 
and user-defined attributes are not limited to C++ implementation. However, the 
FNM will preferably use the same object-oriented programming language used in the 
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reservoir simulation program to facilitate compatibility with the other components of 

the reservoir simulator. 

Novel Use of a OOP Class Hierarchy 

For the computer system using FMN, a novel approach has been taken to 

5 separate much of the member variables from the objects to which they belong. The 
class hierarchy encompasses two sub-hierarchies of classes, one sub-hierarchy 
containing classes to define a variety of "facility types" and the other sub-hierarchy 
containing classes to define a variety of "attribute values." The two sub-hierarchies 
reference each other through their respective base classes - the base facility type class 

10 (FacBase, 200 in Fig. 2) has an association or reference to the base named attribute 
value class (Value, 207 in Fig. 2) via the ValueUse class (206, Fig. 2). In the present 
usage, a facility object will reference many different attribute value objects. Each 
attribute value object contains a specialized unit of data that helps to define the 
characteristics of the facility object. Conceptually, the facility is defined by a 

15 collection of associated objects - the facility object itself, plus all of its associated 
attribute value objects. This class architecture does not strictly follow the 
conventional OOP practice of encapsulation, because a facility's data is not contained 
fully within a single object. However, the attribute value objects are effectively 
encapsulated because their access is still provided via methods of the facility object. 

20 A further unique characteristic of a computer system having a FNM architecture 

is that both the facility type classes and attribute value classes are "generic." 
Specifically, this means that the most specialized classes in each of the two previously 
described sub-hierarchies are capable of representing many different, more 
specialized facility types, and many different attribute values. For example, one type 

25 of specialized facility class is Node (203, Fig. 2). This class is generic because Node 
objects can be instantiated to represent a variety of different node facility types, even 
though the class hierarchy does not have a distinct specialized node class for each 
node facility type. Likewise, the SystemAttributeValue (212, Fig. 2) is a generic 
attribute class because SystemAttributeValue objects can be instantiated to represent a 

30 variety of different attribute types, such as floating point scalars, integer scalars, 
strings, enumerated types, floating point arrays, and integer arrays. 
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A key advantage of the class hierarchy containing generic facility type and 
attribute value classes is that the collection of facility types and attribute values 
supported in the software system can be modified (additional facility types and 
attribute values defined, or existing ones removed) with no change to the class 

5 hierarchy, thereby avoiding the need to recompile source code to accommodate 

changes to facility types or attribute values. In the current system, clear-text file, for 
example, an ASCII file, which can be referred to as a " Data Definitions File" contains 
text entries that define each facility type and all of the attribute values for each facility 
type. This Data Definitions File is processed by a program to create definitions 

10 within the software system. A facility type definition is an object of the FacilityType 
class (302, Fig. 3). An attribute value definition is an object of the 
AttributeDefinition class (305, Fig. 3). When a user works with the system to create a 
simulation model, the definition objects are used by the software whenever the user 
needs to create an instance of a facility, and its accompanying attribute value 

15 instances. 

The FNM overcomes three challenging problems faced by developers of complex 
simulation models of subterranean hydrocarbon-bearing reservoirs and associated 
production facilities: 

1 . Integration of facility and reservoir fluid-flow calculations, which in the past 
20 required models for each that were similar enough to allow for a unified solution 

approach. 

2. The facility fluid-flow calculation model is flexible and easily extensible so that 
the behavior of existing facility types can be extended or additional facility types can 
be added with relatively minor changes to the software. 

25 3. The facility model architecture supports a system of custom facility management 
logic, which is desirable to dynamically control the behavior of facilities during a 
simulation run. 

Facility Management Loeic 

Facility Management Logic ("FML") is a programming system designed for the 
30 operational control of a physical system such as wells and associated surface facilities 
within a reservoir simulation system. The FML comprises three elements: (1) a logic 
interface, referred to herein as "Logic Diagram ", (2) a generator of C++ code (an 
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industry standard programming language) from the Logic Diagram, and (3) a specific 
software design to integrate the FML code with the reservoir simulator itself. 
Logic Diagram 

The Logic Diagram comprises a graphical user interface that a simulator user can 
use to dynamically construct a flow chart of the facility management logic, thereby 
enabling the user to customize simulation of the physical system being simulated. 
The ability to dynamically construct logic refers to the user's ability to construct logic 
as needed to solve facility management problems. Logic can be constructed by the 
user with the logic interface. The logic interface can be used to perform at least one 
step chosen from: developing new logic; selecting and using existing logic that either 
the user or another person previously constructed, perhaps for another facility 
management problem; selecting and modifying existing logic. The flow chart 
comprises icons that represent basic logic control constructs and arrows that indicate 
the progression of logic from one icon to another. The starting point for a Logic 
Diagram consists of a Start icon (representing the entry point to logic), a Finish icon 
(representing the exit point from logic) and a single arrow between the two. Arrows 
have a single direction and they represent the sequence that the icons are evaluated. 
The user's task is to think through the desired logic, and to model that logic by 
building up a representative flow chart. All the logic can be in one high-level 
diagram, or it can be broken up into multiple diagrams, each representing a module of 
logic. Each module of logic is referred to as an Operation. The MainOperation is the 
highest-level Operation. A simplified example of such logic is shown in Fig. 5. 

Fig. 5 diagrammatically illustrates one example of a simplified Logic Diagram 
that shows a simple facility management logic. This logic uses three icon types (three 
Process icons 512, 520 and 526, two Logical Test icons 516 and 522, and one Loop 
icon 514). During the course of performing reservoir simulation, this logic is 
repeatedly executed at regular intervals, each time performing checks and controls on 
producing wells. Specifically, there is a loop over producing wells (Loop icon 514) 
which iterates over all members of the set "producers." In the body of the loop (icons 
516, 520, 522, 526, 528, 530 and the arrows that connect these icons to each other and 
to Loop icon 514) the water rate of a specific producing well in set "producers" is 
checked twice (at Logical Test icons 516 and 522). If the water rate is over 4000 
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barrels per day at Logical Test icon 516, then the well is worked-over (Process icon 
520) in an attempt to reduce the water production. If the water rate is over 5000 
barrels per day at Logical Test icon 522, then the well is closed (Process icon 526). 
Three types of icons form the building blocks of the flow chart. A Process icon is 

5 used to represent equations or calls to operations that act upon facilities to change 
their state in some way (for instance to open or shut-in a well). A Process icon 
represents statement execution, and always has a single inward directing arrow and a 
single outward directing arrow, indicating that there is only one path of execution 
through the Process icon. A Logical Test icon is used to for decision-making in the 

10 flow chart. A Logical Test icon contains a user-specified expression that evaluates to 
true (yes) or false (no). A Logical Test icon has two out-going arrows, one to be 
followed if the expression is true (in Fig. 5, arrow 518 for Logical Test icon 516 and 
arrow 524 for Logical Test icon 522), the other if the expression is false (in Fig. 5, 
arrow 519 for Logical Test icon 516 and arrow 525 for Logical Test icon 522). A 

15 "logical test construct" comprises the Logical Test "icon", and all logic along the Yes 
and No logic branches up to the end marker of the logical test construct (icon 528 for 
Logical Test icon 522 and icon 530 for Logical Test icon 516). A Loop icon is used 
for looping over logic, for example whenever logic must be repeatedly evaluated 
within the flow chart. The Loop icon contains the logic that controls what type of 

20 looping will occur (in Fig. 5, Loop icon 514 is a loop over the members of a set). A 
"loop construct" is comprised of the Loop icon and all logic icons between the Loop 
icon's begin arrow (5 1 5) and next arrow (53 1). Often the facility management logic is 
iterative in nature, and the Loop icon and loop construct are the means to support that 
concept. 

25 The FML user builds up a flow chart for a desired application in a simulation 

model by pointing and clicking in a suitable graphical user interface (a Logic 
Diagram) using a mouse cursor. Using the flow chart displayed in the GUI, the user 
decides which of the three icon types is needed next, and then activates one of three 
toolbar buttons by clicking on the desired button. Next, the user moves the cursor 

30 over an existing arrow in the flow chart that represents the insertion point for the next 
icon. When the user selects the arrow by clicking on it, the system inserts the icon 
and manages the creation and attachment of existing and new arrows. 
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The simulator user's next task is to populate the inserted icon with information 
about specifically what the icon should do. For example, if it is a Logical Test icon, 
then the user defines the expression that is the basis of the true/false test. If it is a 
Loop icon, then information is needed to specify the basis of the iteration, and how 

5 many times the iteration should be done. For a Process icon, the user may need to 
write an equation or specify the details about how an operation is to be invoked. For 
each of these tasks, the user can either type information or can make use of a series of 
context-sensitive menus and dialogs that guide the user through many of the steps 
needed to build-up the contents of an icon. 

10 A specific syntax is preferred for the way information is specified in each icon. 

The syntax can be defined by an FML-specific format referred to as Facility 
Management Control Language. Much of the syntax details are handled 
automatically when the menus and dialogs are used. If the user is typing the 
information directly for each icon, then it is the user's responsibility to adhere to the 

15 syntax rules. If the user prefers not to work with the graphical Logic Diagram, he or 
she can alternately work exclusively in a text-based logic code interface (referred to 
as the Text Mode View) where the user can enter all the Facility Management Control 
Language directly as text, without the visual aides offered within the FML Logic 
Diagram view (icons, arrows, menus or dialogs). Once a user begins to work in the 

20 Text Mode View, he or she cannot return to the Logic Diagram corresponding to the 
logic they are specifying. 

Appendix 1 at the end of this description contains the text from a Text Mode view 
of the logic depicted in flow chart form in Fig. 5. The text appearing in Appendix 1 
conforms to the Facility Management Control Language syntax. The FML 

25 automatically generates this Text Mode view once the user has input the logic using 
the flow chart view. The user does not need to see or work with this view if he or she 
does not want to. On the other hand, some users may choose to enter all the logic 
directly in a Text Mode View, or they may prefer to edit, including entering, 
modifying or deleting the logic once it is in this form. If they do, then the Logic 

30 Diagram view is no longer accessible for this set of logic. As long as the user creates 
or edits logic using only the Logic Diagram view, the user can switch back and forth 
between the two views. 
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Generation of C++ Code 

At the point when the logic is complete, and the user is ready to begin the 
simulation, the system must turn Facility Management Control Language into C++ 
code so that it can be executed with the simulator. If the Logic Diagram was used, 
5 then the system first turns the flow chart into Facility Management Control Language. 
Conversion of logic to object-oriented code may be performed using any suitable 
conversion means. This is an automatic step hidden from the user. If the Text Mode 
was used, then the Facility Management Control Language already in the view is 
used. 

10 The system then begins a process to parse the Facility Management Control 

Language to validate the syntax. As the parsing proceeds, syntactically correct and 
complete pieces of Facility Management Control Language representing the logic is 
automatically converted into ANSI standard object-oriented C++ code. If syntax 
errors are found, they are reported. If there are no errors, then the complete set of 

15 C++ code is written to a file along with the other simulator input data. 

Appendices 2 through 5, located at the end of this description, depict the C++ 
code generated for the sample logic in Fig. 5 and Appendix 1. 

In Appendix 2, the global pointers referring to each facility in the network are 
declared and initialized. This list varies based on the number and type of facilities in 

20 the network, and the facility names specified by the user. The pointers are to the 
interface classes described in later sections, and in Figs. 6 and 7. 

In Appendix 3, all the global facility pointers shown in Appendix 2 are initialized 
to point to the actual facility objects in the simulator. 

In Appendix 4, at each timestep, the simulator invokes the function named 

25 wm_initialize(). This ensures that the global facility pointers refer to the simulator 
objects (by calling wm_initialize_globals() in Appendix 3), and then it calls the 
actual facility management logic. The remaining declarations depict one example of 
dynamic specialization, where the SimWmWelllnterface class is specialized as the 
FunctorO class is declared for the sole purpose of defining a method used to form a set 

30 of producing wells. Dynamic specialization is discussed with respect to Fig. 6, 
discussed in more detail later in this specification. 
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Appendix 5 lists the C++ code for the facility management logic depicted in Fig. 
5 and Appendix 1 . 

When the main simulator is spawned, it extracts the C++ code from the input file 
and puts the code into a file for compilation. The main simulator starts the C++ 

5 compiler that is native to the platform where the simulator will run, and the compiler 
then turns the code into an object file. The C++ linker is then started, to create a 
dynamic link library (on NT) or a shared library (on UNIX). The simulator then loads 
the library, in preparation to execute the code. At the point the main simulator has 
loaded the library, the integrated simulation system is created, comprising the main 

10 simulator and the custom logic. 

At the appropriate point while the integrated simulation system is executing, the 
simulator calls the C++ code corresponding to the MainOperation, which, as 
described earlier, either contains all the facility management logic, or calls to other 
Operations that contain modules of logic. In either case, all the logic is executed, and 

15 upon completion, the code returns control back to the simulator. The simulator then 
continues with its calculations, which are likely impacted by the way in which the 
logic may have changed the state of the facilities. The simulator then continues to 
"move forward" in time to estimate the behavior of the reservoir. Regularly 
throughout this progression, the same facility management logic is executed, thereby 

20 continuously effecting the behavior of the field throughout its simulated life. 
Simulator users can then review the results, and potentially refine the facility 
management logic or the simulator input data. Ultimately, when the results are 
acceptable, the facility management logic can form the basis of the actual facility 
management strategy than will be used in the actual reservoir. 
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Integration with the Reservoir Simulator 

In order for the logic to impact the state of the facilities used by the simulator, it is 
not enough that the simulator execute the logic. The logic, or more specifically the 
object oriented C++ code generated from the logic, is integrated with the main 
simulator's data model. The combination of the main simulation system and the object 
oriented C++ code which is integrated with the main simulator's data model is 
collectively known as the integrated simulation system. In order to integrate the 
object oriented C++ code with the data model, and to make the process more efficient, 
use is made of object-oriented data modeling techniques, specifically multiple 
inheritance and polymorphism, and dynamic specialization. The preferred way in 
which the object-oriented techniques are used will now be described. 
Multiple Inheritance and Polymorphism 

In order to compile and link the C++ code and integrate it with the simulator, the 
C++ code must be aware of various aspects of the simulator's data and methods. This 
is accomplished by "including" files that contain descriptions of the data and methods. 
These files are referred to as header files. These header files are preferably delivered 
with the simulator, so that the compiling and linking can be done. As more header 
files are included in the C++ code, the process of compiling and linking is slowed 
down. To improve performance, it is preferable to reduce the number of header files 
and the data they contain. Multiple inheritance and polymorphism is preferably used 
to accomplish this goal. 

A series of empty "interface" objects are created that represent the simulator 
objects that actually contain the data and methods. This approach is outlined in Fig. 
6 

Fig. 6 schematically illustrates a small portion of the simulator data model. The 
dashed diagonal line 608 divides "Interface" side classes from "Simulator" side 
classes. Solid lines with open triangles at one end indicate inheritance relationships, 
with the base class at the triangle end of the relationship and the derived class at the 
opposite end. The interface classes 600 and 602 on the left of the dashed line 608 are 
essentially empty. They contain no member variables and only empty member 
function declarations ("pure virtual" operations in a C++ implementation of the data 
model). The interface class Welllnterface 602 contains the empty operations 
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"shut_in" and "open". The simulator-side version of Well (606) has the member 
variables for a Well object and the real implementation for the member functions 
"shut_in" and "open". Multiple inheritance permits the use and declaration of the 
interface classes. The simulator- side Well class 606 inherits both from the interface 
class Welllnterface 602 via inheritance arrow 609 and from the simulator class Object 
604 via inheritance arrow 610. Welllnterface class 602 is a specialization of 
Objectlnterface 600. Simulator-side base class Object 604 inherits from 
Objectlnterface class 600 via inheritance arrow 611. Polymorphism permits the 
empty interface operation "open" to be called on the interface class Welllnterface 
602, and then to execute the real "open" operation on the actual simulator class Well 
606. This same paradigm is used to represent every simulator facility class with a 
corresponding interface class. 

Only the interface objects are needed to be able to compile and link the logic. 
Since the compiler sees only the simple interface classes, complexity and proprietary 
information remains hidden from the compiler and from anyone who may be looking 
at the files. The Facility Management Control Language shown in Appendix 2 
contains only interface object pointers, not simulator object pointers. Since the 
interface objects are empty, the compile/link process is fast, and details about the set 
of data and methods needed by the simulator are not disclosed (for example to anyone 
who may be scanning the computer disk for files). Multiple inheritance allows the 
interface classes to be related to the simulator classes, and polymorphism allows real 
methods to be invoked indirectly using empty methods declared on the interface 
classes. 

Dynamic Specialization 

Each logic Operation developed by a simulator user is associated either with the 
facility network in general or with a specific facility. The choice is made by the user 
depending on the nature of the logic. In order for these Operations to be executable 
within the simulator data model framework already described, the Operations need to 
become methods of classes. Since C++ does not permit adding methods to existing 
classes without recompiling the simulator (or significant portions thereof), the 
methods have to be added to new classes that only the logic code is aware of. This 
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leads to the design wherein at runtime new classes are specialized from the existing 
interface classes. This approach is schematically shown in Fig. 7. 

Fig. 7, which extends the diagram in Fig. 6, illustrates how dynamic specialization 
is used to support user-written operations. The diagonal dashed line 710 separates the 

5 interface classes 700, 702 and 704 from the simulator classes 706 and 708. Here, the 
user has written a new operation named "special_well_logic" for facilities of type 
Well. To support the declaration of new operations, the interface classes are 
automatically specialized at runtime using class declarations that are written by the 
Facility Management Control Language-to-C++ parser (hence the term "dynamic 

10 specialization"). In this example, the class "WellInterface_User" 704 is created by the 
parser, and it contains the declaration of and the code comprising the new logic 
operation. Therefore, both the built-in operations (from Welllnterface 702 and 
Objectlnterface 700) and user operations (from WellInterface_User 704) can be called 
for facilities that are of type Well. 

15 The data model is extended, but only for the logic code. The rest of the simulator 

is unaware of the extension, so from the simulator's perspective the data model has 
not changed, and therefore the simulator code does not have to be compiled in order 
to use the user-defined logic. Meanwhile, the logic is aware of the new classes and 
recognizes the new operations. From the user's perspective, built-in operations that 

20 the simulator is aware of, and custom operations that the simulator is unaware of, are 
invoked exactly the same way. Furthermore, the new operations, which belong to 
classes that are derived from the interface classes, have access to the same built-in 
operations and data that the interface classes do. The result is an apparent seamless 
integration between simulator code and custom code. 

25 The FML system can be integrated by persons skilled in the art with a reservoir 

simulator that uses object-oriented design, preferably coded using C++. The 
prerequisite for the FML system is that it must be integrated with a simulator that has 
an object-oriented design. 

It is the object-oriented principles that support the concepts of multiple 

30 inheritance, polymorphism and class specialization, which are necessary to seamlessly 
integrate FML with a simulator. The FML can be used with a simulator that is 
implemented using an object-oriented language other than C++. Other object- 
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oriented programming ("OOP") techniques could be used, including without 
limitation, Simula, Eiffel, Smalltalk, Objective-C, or variations thereof The practice 
of FML is not limited to a particular OOP language. 

If another OOP system other than C++ is used, the C++ code generation step 
5 would need to be revised to generate the alternative OOP language, and this also 
implies that the simulator would also have an object oriented data model. 

The reservoir simulation model of this invention allows the user to develop 
custom facility management logic to simulate practically any strategy. As such, the 
permutations of logic are practically infinite. 
10 Flexible and Extensible Data Model Design 

The reservoir simulation system preferably also uses object-oriented features of 
C++ to address the problem of extensibility. The software design preferably uses a 
collection of C++ classes (referred to collectively as a "data model") to define objects 
with prescribed data (member variables) and behavior (member functions). The 
15 member variables and member functions are hard-coded, that is, they are defined 
within the classes themselves, meaning that making extensions to the classes; for 
example, new member variables, new member functions, new classes, requires code 
changes and recompilation. 

In a typical OOP class hierarchy, all member variables for a given class are 
20 encapsulated within that class. Therefore, objects that are instances of that class 

protect and control access to their member variables through their public interfaces, 
which are the member functions. 

The principle of the new reservoir simulation system and the best mode 
contemplated for applying that principle have been described. It will be apparent to 
25 those skilled in the art that various changes may be made to the embodiments 
described above without departing from the spirit and scope of this reservoir 
simulation system as described above. It is, therefore, to be understood that this 
invention is not limited to the specific details shown and described. 
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Appendix 1 

// Start of Global Variable Declarations 
// End of Global Variable Declarations 
5 wm_main ( ) 

{ // Start of module code 

// Start of Local Variable Declarations 
Set<Well> producers ; 
Ref<Well> well_in_set ; 
10 // End of Local Variable Declarations 

{ // Start of Process icon construct 

producers. form_set (candidate. facility_descriptor == "producer") ; 
} // End of Process icon construct 
{ // Start of loop construct 
15 foreach well_in_set in producers 

{ // Start of loop body 

{ // Start of logical test construct 

if ( ( well_in_set.get_phase_rate (WATER) > 4000. ) ) 

{ 

20 { // Start of Process icon construct 

well_in_set . workover (WRKOVR_CLOSE, BY_RESCONN, 
"WCUT>.6"); 
} // End of Process icon construct 
{ // Start of logical test construct 
25 if ( ( well_in_set.get_phase_rate (WATER) > 5 000. ) ) 

{ 

{ // Start of Process icon construct 

well_in_set . shut_in (NO_CROSS_FLOW) ; 
) // End of Process icon construct 

30 } 

else 

{ 
} 

} // End of logical test construct 

35 } 

else 

{ 

} 

} // End of logical test construct 
40 } // End of loop body 

} // End of loop construct 

{ 

return ; 
} 
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Appendix 2 



♦include "q: \rell_sim\bin\runtime_includes\wm_headers.h" 

SimWmWellNodelnterf ace* wellnode_potnorth = NULL; 
SimWmWellNode Interface* wellnode_potsouth = NULL; 
SimWmWellNodelnterf ace* wellnode_prodlOO = NULL; 
SimWmWellNodelnterf ace* wellnode_prod82 = NULL; 
SimWmNetworkNodelnterf ace* netnode_inj ector = NULL; 
SimWmNetworkNodelnterface* netnode_producer = NULL; 
SimWmTerminalNodelnterf ace* termnode_inj ector = NULL; 
SimWmTerminalNodelnterf ace* termnode_producer = NULL; 
SimWmReservoirConnectionlnterf ace* potnorth_rc_0 = NULL; 
SimWmReservoirConnectionlnterf ace* potnorth_rc_l = NULL; 
SimWmReservoirConnectionlnterf ace* potnorth_rc_2 = NULL; 
SimWmReservoirConnectionlnterf ace* potsouth_rc_0 = NULL; 
SimWmReservoirConnectionlnterface* potsouth_rc_l = NULL; 
SimWmReservoirConnectionlnterface* potsouth_rc_2 = NULL; 
SimWmReservoirConnectionlnterf ace* prodlOO_rc_0 = NULL; 
SimWmReservoirConnectionlnterf ace* prodlOO_rc_l = NULL; 
SimWmReservoirConnectionlnterf ace* prodl00_rc_2 = NULL; 
SimWmReservoirConnectionlnterf ace* prod82_rc_0 = NULL; 
SimWmReservoirConnectionlnterf ace* prod82_rc_l = NULL; 
SimWmReservoirConnectionlnterf ace* prod82_rc_2 = NULL; 
SimWmTerminalConnectionlnterf ace* termconn_inj ector = NULL 
SimWmNetworkConnectionlnterf ace* netconn_potnorth = NULL; 
SimWmNetworkConnectionlnterf ace* netconn_pot south = NULL; 
SimWmNetworkConnectionlnterf ace* netconn_j>rodlOO = NULL; 
SimWmTerminalConnectionlnterface* termconn_producer = NULL 
SimWmNetworkConnectionlnterf ace* netconn_prod82 = NULL; 
SimWmWelllnterface* potnorth = NULL; 
SimWmWellBranchlnterface* wellbranchl = NULL; 
SimWmWelllnterface* potsouth = NULL; 
SimWmWellBranchlnterface* wellbranch2 = NULL; 
SimWmWelllnterface* prodlOO = NULL; 
SimWmWellBranchlnterface* wellbranch3 = NULL; 
SimWmWelllnterface* prod82 = NULL; 
SimWmWellBranchlnterface* wellbranch4 = NULL; 
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Appendix 3 



void wm_initialize_globals ( ) 

{ 

static int wmig_was_called = 0; 
5 if (wmig_was_called) 

return; 
wmig_was_called = 1; 

wellnode_potnorth = (SimWmWellNodelnterf ace* ) 

SimWmFNControllerUtility: : lookup_obj ect ( "wellnode_potnorth" , 
10 VM_GetInterfaceID<SimWmWellNodeInterface>() () ) ; 

wellnode_potsouth = (SimWmWellNodelnterf ace* ) 

SimWmFNControllerUtility: : lookup_object ( "wellnode_potsouth" , 
VM_GetInterfaceID<SimWmWellNodeInterface>() ( ) ) ; 
wellnode_prodl00 = (SimWmWellNodelnterf ace*) 
15 SimWmFNControllerUtility: : lookup_object ("wellnode_prodl00" / 

VM_GetInterfaceID<SimWmWellNodeInterface>() ( ) ) ; 
wellnode_prod82 = (SimWmWellNodelnterf ace* ) 

SimWmFNControllerUtility: : lookup_ob j ect ( "wellnode_prod82 " , 
VM_GetInterfaceID<SimWmWellNodeInterface>() ( ) ) ; 
20 netnode_inj ector = ( SimWmNetworkNodelnterf ace* ) 

SimWmFNControllerUtility: : lookup_object ("netnode_inj ector", 
VM_GetInterfaceID<SimWmNetworkNodeInterface> () () ); 
netnode_producer = (SimWmNetworkNodelnterf ace* ) 

SimWmFNControllerUtility: : lookup_ob j ect ( "netnode_producer " , 
25 VM_GetInterfaceID<SimWmNetworkNodeInterface> () ()); 

termnode_inj ector = (SimWmTerminalNodelnterf ace* ) 

SimWmFNControllerUtility: : lookup_ob j ect ( "termnode_inj ector" , 
VM_GetInterfaceID<SimWmTerminalNodeInterface> () ( ) ) ; 
termnode_producer = (SimWmTerminalNodelnterf ace*) 
30 SimWmFNControllerUtility: : lookup_obj ect ( "termnode_producer " , 

VM_GetInterf aceID<SimWmTerminalNodeInterf ace> ( ) ( ) ) ; 
potnorth_rc_0 = (SimWmReservoirConnectionlnterf ace*) 

SimWmFNControllerUtility: : lookup_obj ect ( "potnorth_rc_0" , 
VM_GetInterf aceID<SimWmReservoirConnectionInterf ace> ( ) ( ) ) ; 
35 potnorth_rc_l = (SimWmReservoirConnectionlnterf ace* ) 

SimWmFNControllerUtility: : lookup_object ( "potnorth_rc_l" , 
VM_GetInterfaceID<SimWmReservoirConnectionInterface>() () ) ; 
potnorth_rc_2 = (SimWmReservoirConnectionlnterf ace* ) 

SimWmFNControllerUtility: : lookup_obj ect ("potnorth_rc_2" , 
40 VM_GetInterf aceID<SimWmReservoirConnectionInterf ace> { ) ()); 

potsouth_rc_0 = (SimWmReservoirConnectionlnterf ace* ) 

SimWmFNControllerUtility: : lookup_object ("potsouth__rc_0", 
VM_GetInterf aceID<SimWmReservoirConnectionInterf ace> ( ) ( ) ) ; 
potsouth_rc_l = (SimWmReservoirConnectionlnterf ace*) 
45 SimWmFNControllerUtility: : lookup_object ("potsouth_rc_l" , 

VM_GetInterf aceID<SimWmReservoirConnectionInterf ace> ( ) ( ) ) ; 
potsouth_rc_2 = (SimWmReservoirConnectionlnterf ace* ) 

SimWmFNControllerUtility: : lookup_obj ect ("potsouth_rc_2" , 
VM_GetInterf aceID<SimWmReservoirConnectionInterf ace> ( ) ( ) ) ; 
50 prodl00_rc_0 = (SimWmReservoirConnectionlnterf ace*) 

SimWmFNControllerUtility: : lookup_object ("prodl00_rc_0", 

VM_GetInterfaceID<SimWmReservoirConnectionInterf ace> () () ) 
prodl00_rc_l = (SimWmReservoirConnectionlnterf ace*) 

SimWmFNControllerUtility: : lookup_obj ect ("prodl00_rc_l" , 
55 VM_GetInterfaceID<SimWmReservoirConnectionInterface> ( ) ())/ 

prodl00_rc_2 = (SimWmReservoirConnectionlnterf ace*) 

SimWmFNControllerUtility: : lookup_ob j ect ( "prodl00_rc_2 " , 
VM_GetInterf aceID<SimWmReservoirConnectionInterf ace> ( ) ( ) ) ; 
prod82_rc_0 = (SimWmReservoirConnectionlnterf ace* ) 
60 SimWmFNControllerUtility: : lookup_obj ect ("prod82_rc_0", 

VM_GetInterf aceID<SimWmReservoirConnectionInterf ace> ( ) () ) ; 
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prod82_rc_l = (SimWmReservoirConnectionlnterf ace*) 

SimWmFNControllerUtility : : lookup_ob j ect ("prod82_rc_l", 
VM_GetInterfaceID<SimWniReservoirConnectionInterface>() () ) ; 
prod82_rc_2 = (SimWmReservoirConnectionlnterface*) 
5 SimWmFNControllerUtility: : lookup_ob j ect ( "prod8 2_rc_2 " , 

VM_GetInterf aceID<SimWmReservoirConnectionInterf ace> ( ) () ) ; 
termconn_injector = (SimWmTerminalConnectionlnterf ace* ) 

SimWmFNControllerUtility: : lookup_ob j ect ( "termconn_inj ector 
VM_GetInterf aceID<SimWmTerminalConnectionInterface> ( ) ()) ; 
10 netconn_potnorth = (SimWmNetworkConnectionlnterf ace* ) 

SimWmFNControllerUtility: : lookup_obj ect ( "netconn_potnorth" 
VM_GetInterf aceID<SimWmNetworkConnectionInterf ace> () () ) ; 
netconn_potsouth = (SimWmNetworkConnectionlnterf ace*) 

SimWmFNControllerUtility: : lookup_ob j ect ( "net conn_pot south" 
15 VM_GetInterfaceID<SimWmNetworkConnectionInterface> () ()); 

netconn_prodl00 = (SimWmNetworkConnectionlnterf ace*) 

SimWmFNControllerUtility: : lookup_obj ect ( "netconnjprodlOO " , 
VM_GetInterfaceID<SimWmNetworkConnectionInterf ace> () () ) ; 
termconn_producer = (SimWmTerminalConnectionlnterf ace* ) 
20 SimWmFNControllerUtility: : lookup_obj ect ( "termconn_producer 

VM_GetInterf aceID<SimWmTerminalConnectionInterf ace> ( ) () ) ; 
netconn_prod82 = (SimWmNetworkConnectionlnterf ace* ) 

SimWmFNControllerUtility: : lookup_ob j ect ( "netconn_prod82" , 
VM_GetInterf aceID<SimWmNetworkConnectionInterf ace> ( ) () ) ; 
25 potnorth = (SimWmWelllnterf ace*) 

SimWmFNControllerUtility: : lookup_obj ect ("potnorth", 
VM_GetInterfaceID<SimWmWellInterface> () ( ) ) ; 
wellbranchl = (SimWmWellBranchlnterf ace* ) 

SimWmFNControllerUtility: :lookup_obj ect ("wellbranchl" , 
30 VM_GetInterf aceID<SimWmWellBranchInterf ace> ( ) ()); 

potsouth = (SimWmWelllnterf ace*) 

SimWmFNControllerUtility: : lookup_ob j ect ("potsouth", 
VM_GetInterfaceID<SimWmWellInterface>() () ) ; 
wellbranch2 = (SimWmWellBranchlnterf ace* ) 
35 SimWmFNControllerUtility: : lookup_ob j ect ("wellbranch2", 

VM_GetInterf aceID<SimWmWellBranchInterf ace> ( ) () ) ; 
prodlOO = (SimWmWelllnterf ace*) 

SimWmFNControllerUtility: : lookup_obj ect ("prodlOO", 
VM_GetInterfaceID<SimWmWellInterface>() ()) ; 
40 wellbranch3 = (SimWmWellBranchlnterf ace* ) 

SimWmFNControllerUtility: : lookup_object ("wellbranch3", 
VM_GetInterf aceID<SimWmWellBranchInterface> ( ) ( ) ) ; 
prod82 = (SimWmWelllnterf ace*) 

SimWmFNControllerUtility: : lookup_obj ect ("prod82", 
45 VM_GetInterf aceID<SimWmWellInterf ace> ( ) ( ) ) ; 

wellbranch4 = (SimWmWellBranchlnterf ace* ) 

SimWmFNControllerUtility: : lookup_ob j ect ( "wellbranch4 ", 
VM_GetInterfaceID<SimWmWellBranchInterface> ( ) () ) ; 
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Appendix 4 

void wm_main() ; 
extern "C" { 
5 void wm_initialize () 

{ 

wm_initialize_globals ( ) ; 
wm_main ( ) ; 

} 

10 } 

class FunctorO : public Selector<SimWmWellInterf ace> 
{ 

public : 

FunctorO ( ) ; 
15 virtual -FunctorO (); 

virtual Bool operator ( ) (SimWmWelllnterf ace* candidate) const; 

>; 

FunctorO : : FunctorO ( ) 

{ 

20 //do nothing else 

} 

FunctorO : : ~FunctorO ( ) 
{ 

//do nothing else 

25 } 

Bool FunctorO: : operator () (SimWmWelllnterf ace* candidate) const 
{ 

return (*( (String*) (candidate- 
>lookup_attribute ( "f acility_descriptor" , 0, 
30 "~ TYPE_String) ) ) == "producer"); 

} 



EXPRESS MAIL NO. ET607889812US 
DECEMBER 6, 2001 



32 



Appendix 5 



void wm_main ( ) 
5 { 

Set<SimWmWellInterface> producers; 
SimWmWelllnterface* well_in_set = NULL; 

{ 

producers . f orm_set (new FunctorO ( ) ) ; 

10 } 

{ 

producers . reset ( ) ; 

while (producers . get_next (well_in_set) ) 

{ 

15 PhaseTypeEnum scratch 0 = (PhaseTypeEnum) (WATER) ; 

if ( (well_in_set->get_phase_rate (scratch 0) > 4000.)) 

{ 

{ 

WorkoverTypeEnum scratch 1 = 

20 (WorkoverTypeEnum) (WRKOVR_CLOSE) ; 

WorkoverLevelTypeEnum scratc h 2 = 

(WorkoverLevelTypeEnum) (BY_RESCONN) / 
String scratch 3 = "WCUT>.6"; 

well_in_set->workover (scratch 1, scratch 2 , 



25 scratch 3) ; 



30 5000.)) 



} 



PhaseTypeEnum scratch 4 = (PhaseTypeEnum) (WATER) ; 

if ( (well_in_set->get_phase_rate (scratch 4) > 



ShutlnTypeEnum scratch 5 = 

ShutlnTypeEnum) (NO_CROSS_FLOW) ; 



. set->shut_in (scratch 5); 



else 

{ 
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Appendix 6 

This excerpt from the Data Definitions File illustrates the hierarchical nature of 
facility type definitions and named attribute definitions. The facility types included in 
this excerpt demonstrate the base facility type levels needed to define the end-user 
facility type "NetworkNode." The excerpt includes the definitions for three named 
attributes that a NetworkNode would have. Two of the named attributes are defined at 
the FacObject level and the third is defined at the Node level. 

// Comments in the Data Definitions File begin with "//" 

// "\" at the end of a line indicates continuation on the subsequent line 
// Highest level base facility type is FacBase 



FacBase" 

"FacBase Definition" 



"FacBase" 
"View Name" 



"restart rule" 
"simulator rule" 
"user rule" 
"valid rule" 



Definition: 
Is Obsolete: 
IOArchive Class Name: 
Long Name: 
Vi ew Name : 
Decomp Hierarchy: 
Dimensionality: 
Req by Restart Rule: 
Req by Simulator Rule 
Req by User Rule: 
Valid Rule: 
Bitmap Resource ID: 0 
Instantiable: F 
Member of Granule: F 
Simulator Variable: F 
Time Variant: F 
Decomposable: F 
Spec for all Nodes: F 
Change Invalidates : F 
Parent Facility: 
Is Abstract: T 

Parent type of FacObject is FacBase 
s: "FacObject" 

Definition: 

Is Obsolete: 

IOArchive Class Name: 

Long Name: 

View Name: 

Decomp Hierarchy: 

Dimensionality: 

Req by Restart Rule: 

Req by Simulator Rule: 

Req by User Rule: 

Valid Rule: 

Bitmap Resource ID: 

Instantiable: 

Member of Granule: 

Simulator Variable: 

Time Variant: 

Decomposable : 

Spec for all Nodes: 

Change Invalidates: 

Parent Facility: 

Is Abstract: 



'FacObject Definiti 



"FacObject" 
"View Name" 



'restart rule" 
'simulator rule" 
'user rule" 
'valid rule" 



"FacBase" 



// is_active is a named attribute defined at the FacObject level 



Attribute Name: 



"is active" 
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Description: 

Unit Group Name: 

Default Value: 

Enmti Type: 

Is Array: 

Is Fixed Length: 

Number of Values: 

Dimensionality: 

Attribute Type: 

Is Class Attribute: 

Is Visible to User: 

Is Time Variant: 

Is Results Only: 

Is Used by Simulator: 



"is_active" 



// f acility_descriptor : 
FacObject level 

Attribute Name: 

Description: 



be i 



sed" 



. named attribute defined at the 



ty_descriptor" 

"User-defined string that < 



to logically group this 
with similar facilities" 



Unit Group Name: 

Default Value: 

Enum Type: 

Is Array: 

Is Fixed Length: 

Number of Values : 

Dimensionality: 

Attribute Type: 

Is Class Attribute: 

Is Visible to User: 

Is Time Variant: 

Is Results Only: 

Is Used by Simulator: 



FacObj ect 
Node " 

"Node Definitio 

F 

"NA" 

"View Name" 



"restart rule" 
"simulator rule" 
"user rule" 
"valid rule" 



// Parent type of Node 
Name : 

Definition: 
Is Obsolete: 
IOArchive Class Name: 
Long Name : 
View Name: 
Decomp Hierarchy: 
Dimensionality: 
Reg by Restart Rule: 
Req by Simulator Rule 
Req by User Rule: 
Valid Rule: 
Bitmap Resource ID: 0 
Instantiable: F 
Member of Granule: F 
Simulator Variable: F 
Time Variant: F 
Decomposable: F 
Spec for all Nodes: F 
Change Invalidates: F 
Parent Facility: "FacObject" 
Is Abstract: T 

// depth_node is a named attribute defined at the Node level 
Attribute Name: "depth_node" 

Description: "Depth of the node" 
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Unit Group Name: 

Default Value: 

Enum Type: 

Is Array: 

Is Fixed Length: 

Number of Values: 

Dimens i onal i ty : 

Attribute Type: 

Is Class Attribute: 

Is Visible to User: 

Is Time Variant: 

Is Results Only: 

Is Used by Simulator: 



"Depths' 
"0.0" 



Note that the FNM C++ 



// Parent type of GeneralFacilityNode is Node, 
data model 

// does not contain this class - it is only defined in the Data 
Definitions File. 

Nam e: "GeneralFacilityNode" 

Definition: "generic class used to represent most 

facility types" 

Is Obsolete: F 

IOArchive Class Name: 

Long Name: 

View Name: 

Decomp Hierarchy: 

Dimensionality: 

Req by Restart Rule: 

Req by Simulator Rule: 

Req by User Rule: 

Valid Rule: 

Bitmap Resource ID: 0 

Instantiable: F 

Member of Granule: F 

Simulator Variable: F 

Time Variant: F 

Decomposable: F 

Spec for all Nodes: F 

Change Invalidates: F 

Parent Facility: "Node" 

Is Abstract: T 

// Parent type of NetworkNode is GeneralFacilityNode. Note that the FNM 
C++ data 

// model does not contain this class - it is only defined in the Data 
Definitions 
// File. 

Name: "NetworkNode" 



"GeneralFacilityNode" 
"View Name" 



"restart rule" 
"simulator rule" 
"user rule" 
"valid rule" 



Definition: 

Is Obsolete: 

IOArchive Class Name: 

Long Name: 

View Name: 

Decomp Hierarchy: 

Dimensionality: 

Req by Restart Rule: 

Req by Simulator Rule: 

Req by User Rule: 

Valid Rule: 

Bitmap Resource ID: 

Instantiable: 

Member of Granule: 

Simulator Variable: 

Time Variant: 



"NetworkNode Definition" 



"restart rule" 
"simulator rule" 
"user rule" 
"valid rule" 
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Decomposable: F 

Spec for all Nodes: F 

Change Invalidates: F 

Parent Facility: "GeneralFacilityNode" 

Is Abstract: F 
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